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1 Introduction 

The theory of polyhedral surfaces and, more generally, the field of discrete differ- 
ential geometry are presently emerging on the border of differential and discrete 
geometry Whereas classical differential geometry investigates smooth geometric 
shapes (such as surfaces), and discrete geometry studies geometric shapes with 
a finite number of elements (polyhedra), the theory of polyhedral surfaces aims 
at a development of discrete equivalents of the geometric notions and methods 
of surface theory. The latter appears then as a limit of the refinement of the 
discretization. Current progress in this field is to a large extent stimulated by 
its relevance for computer graphics and visualization. 

One of the central problems of discrete differential geometry is to find proper 
discrete analogues of special classes of surfaces, such as minimal, constant mean 
curvature, isothermic surfaces, etc. Usually, one can suggest various discretiza- 
tions with the same continuous limit which have quite different geometric prop- 
erties. The goal of discrete differential geometry is to find a discretization which 
inherits as many essential properties of the smooth geometry as possible. 

Our discretizations are based on quadrilateral meshes, i.e. we discretize pa- 
rameterized surfaces. For the discretization of a special class of surfaces, it is 
natural to choose an adapted parameterization. In this paper, we investigate 
conformal discretizations of surfaces, i.e. discretizations in terms of circles and 
spheres, and introduce a new discrete model for minimal surfaces. See Figs. H 
and In comparison with direct methods (see, in particular, |23j). leading 
usually to triangle meshes, the less intuitive discretizations of the present paper 
have essential advantages: they respect conformal properties of surfaces, possess 
a maximum principle (see Remark on p. etc. 

We consider minimal surfaces as a subclass of isothermic surfaces. The anal- 
ogous discrete surfaces, discrete S-isothermic surfaces |1] consist of touching 
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Figure 1: A discrete minimal Ennepcr surface (left) and a discrete minimal 
catenoid (right). 

spheres, and of circles which intersect the spheres orthogonally in their points 
of contact. See Fig. ^ (right). Continuous isothermic surfaces allow a duality 
transformation, the Christoffel transformation. Minimal surfaces are character- 
ized among isothermic surfaces by the property that they are dual to their Gauss 
map. The duality transformation and the characterization of minimal surfaces 
carries over to the discrete domain. Thus, one arrives at the notion of discrete 
minimal S-isothermic surfaces, or discrete minimal surfaces for short. The role 
of the Gauss maps is played by discrete S-isothermic surfaces, the spheres of 
which all intersect one fixed sphere orthogonally. Due to a classical theorem of 
Koebe (see Section El any 3-dimcnsional combinatorial convex polytope can be 
(essentially uniquely) realized as such a Gauss map. 

This definition of discrete minimal surfaces leads to a construction method 
for discrete S-isothermic minimal surfaces from discrete holomorphic data, a 
form of a discrete Weierstrass representation (see Section |5j). Moreover, the 
classical "associated family" of a minimal surface, which is a one-parameter 
family of isometric deformations preserving the Gauss map, carries over to the 
discrete setup (see Section EJl. 

Our general method to construct discrete minimal surfaces is schematically 
shown in the following diagram. (See also Fig|15l1 

continuous minimal surface 
image of curvature lines under Gauss-map 
cell decomposition of (a branched cover of) the sphere 
orthogonal circle pattern 
Koebe polyhedron 

A 

discrete minimal surface 

As usual in the theory on minimal surfaces |18j . one starts constructing such a 
surface with a rough idea of how it should look. To use our method, one should 
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Figure 2: A discrete minimal Schwarz P-surface (left) and a discrete minimal 
Scherk tower (right). 

understand its Gauss map and the combinatorics of the curvature line pattern. 
The image of the curvature line pattern under the Gauss map provides us with a 
cell decomposition of (a part of) S 2 or a covering. From these data, applying the 
Koebe theorem, we obtain a circle packing with the prescribed combinatorics. 
Finally, a simple dualization step yields the desired discrete minimal surface. 

Let us emphasize that our data, besides possibly boundary conditions, are 
purely combinatorial — the combinatorics of the curvature line pattern. All faces 
are quadrilaterals and typical vertices have four edges. There may exist dis- 
tinguished vertices (corresponding to the ends or umbilic points of a minimal 
surface) with a different number of edges. 

The most nontrivial step in the above construction is the third one listed in 
the diagram. It is based on the Koebe theorem. It implies the existence and 
uniqueness for the discrete minimal S-isothcrmic surface under consideration, 
but not only this. This theorem can be made an effective tool in constructing 
these surfaces. For that purpose, we use a variational principle from 0, |28| for 
constructing circle patterns. This principle provides us with a variational de- 
scription of discrete minimal S-isothermic surfaces and makes possible a solution 
of some Plateau problems as well. 

In Section [7J we prove the convergence of discrete minimal S-isothcrmic 
surfaces to smooth minimal surfaces. The proof is based on Schramm's approx- 
imation result for circle patterns with the combinatorics of the square grid |2fij . 
The best known convergence result for circle patterns is C°° -convergence of 
circle packings |14j . It is an interesting question whether the convergence of 
discrete minimal surfaces is as good. 

Because of the convergence, the theory developed in this paper may be used 
to obtain new results in the theory of smooth minimal surfaces. A typical 
problem in the theory of minimal surfaces is to decide whether surfaces with 
some required geometric properties exist, and to construct them. The discovery 
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of the Costa-Hoffman-Meeks surface a turning point of the modern theory 
of minimal surfaces, was based on the Weierstrass representation. This powerful 
method allows the construction of important examples. On the other hand, it 
requires a specific study for each example; and it is difficult to control the 
embeddedness. Kapouleas |23 proved the existence of new embedded examples 
using a new method. He considers finitely many catcnoids with the same axis 
and planes orthogonal to this axis and shows that one can desingularize the 
circles of intersection by deformed Schcrk towers. This existence result is very 
intuitive, but it gives no lower bound for the genus of the surfaces. Although 
some examples with lower genus are known (the Costa-Hoffman-Meeks surface 
and generalizations |2()j). proving the existence of Kapouleas' surfaces with given 
genus and to construct them using conventional methods is very difficult |3U| . 
Our method may become helpful to address these problems. At the present time, 
however, the construction of new minimal surfaces from discrete ones remains 
a challenge. 

Apart from discrete minimal surfaces, there are other interesting subclasses 
of S-isothcrmic surfaces. In future publications, we plan to treat discrete con- 
stant mean curvature surfaces in Euclidean 3-space and Bryant surfaces |7] , |10j . 
(Bryant surfaces are surfaces with constant mean curvature 1 in hyperbolic 3- 
space.) Both are special subclasses of isothermic surfaces that can be charac- 
terized in terms of surface transformations. (See 0] and for a definition 
of discrete constant mean curvature surfaces in 1R 3 in terms of transformations 
of isothermic surfaces. See |17) for the characterization of continuous Bryant 
surfaces in terms of surface transformations.) 

More generally, we believe that the classes of discrete surfaces considered in 
this paper will be helpful in the development of a theory of discrete conformally 
parameterized surfaces. 

2 Discrete S-isothermic surfaces 

Every smooth immersed surface in 3-space admits curvature line parameters 
away from umbilic points, and every smooth immersed surface admits conformal 
parameters. But not every surface admits a curvature line parameterization that 
is at the same time conformal. 

Definition 1. A smooth immersed surface in R 3 is called isothermic if it admits 
a conformal curvature line parameterization in a neighborhood of every non- 
umbilic point. 

Geometrically, this means that the curvature lines divide an isothermic sur- 
face into infinitesimal squares. An isothermic immersion (a surface patch in 
conformal curvature line parameters) 

/ : R 2 D D -> R 3 

(x,y) h-> f(x,y) 

is characterized by the properties 

\\fx\\ = \\fy\\, fx-Lfy, fxy& S P an{ f x , f y } . (1) 

Being an isothermic surface is a Mobius-invariant property: A Mobius trans- 
formation of Euclidean 3-spacc maps isothermic surfaces to isothermic surfaces. 
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The class of isothermic surfaces contains all surfaces of revolution, all quadrics, 
all constant mean curvature surfaces, and, in particular, all minimal surfaces 
(see Theorem^J. In this paper, we are going to find a discrete version of mini- 
mal surfaces by characterizing them as a special subclass of isothermic surfaces 
(see Section 0J. 

While the set of umbilic points of an isothermic surface can in general be 
more complicated, we are only interested in surfaces with isolated umbilic points, 
and also in surfaces all points of which arc umbilic. In the case of isolated umbilic 
points, there are exactly two orthogonally intersecting curvature lines through 
every non-umbilic point. An umbilic point has an even number 2k (k ^ 2) 
of curvature lines originating from it, evenly spaced at n/k angles. Minimal 
surfaces have isolated umbilic points. If, on the other hand, every point of the 
surface is umbilic, then the surface is part of a sphere (or plane) and every 
conformal parameterization is also a curvature line parameterization. 

Definition of discrete isothermic surfaces was already suggested in |2j . 
Roughly speaking, a discrete isothermic surface is a polyhedral surface in 3- 
space all faces of which are conformal squares. To make this more precise, we 
use the notion of a "quad-graph" to describe the combinatorics of a discrete 
isothermic surface, and we define "conformal square" in terms of the cross-ratio 
of four points in R 3 . 

A cell decomposition T> of an oriented two-dimensional manifold (possibly 
with boundary) is called a quad-graph, if all its faces are quadrilaterals, that 
is, if they have four edges. The cross-ratio of four points Z\, z 2 , Z4 in the 
Riemann sphere (D = (DU{oo}is 

/ s (^l-^2)(^3-^4) 

\Z2 - zz){z A - Zl) 

The cross-ratio of four points in 1R 3 can be defined as follows: Let S be a sphere 
(or plane) containing the four points. S is unique except when the four points lie 
on a circle (or line). Choose an orientation on S and an orientation preserving 
conformal map from S to the Riemann sphere. The cross-ratio of the four points 
in R 3 is defined as the cross-ratio of the four images in the Riemann sphere. 
The two orientations on S lead to complex conjugate cross-ratios. Otherwise, 
the cross-ratio does not depend on the choices involved in the definition: neither 
on the conformal map to the Riemann sphere, nor on the choice of S when the 
four points lie in a circle. The cross-ratio of four points in R 3 is thus defined 
up to complex conjugation. (For an equivalent definition involving quaternions, 
see [3], HS|.) The cross-ratio of four points in R 3 is invariant under Mobius 
transformations of R 3 . Conversely, if pi, p 2 , P3, Pa € R 3 have the same cross- 
ratio (up to complex conjugation) as p[, p' 2 , p'3, p\ € R 3 , then there is a Mobius 
transformation of R 3 which maps each pj to p'j . 

Four points in R 3 form a conformal square, if their cross-ratio is —1, that is, 
if they are Mobius-equivalent to a square. The points of a conformal square lie 
on a circle (see Fig- EJ • 

Definition 2. Let T> be a quad-graph such that the degree of every interior 
vertex is even. (That is, every vertex has an even number of edges.) Let V(T>) 
be the set of vertices of V. A function 

/ : V{V) -> R 3 
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Figure 3: Left: A conformal square. The sides a, a', b, b' are interpreted as 
complex numbers. Right: Right-angled kites are conformal squares. 

is called a discrete isothermic surface, if for every face of T> with vertices V\, 
V2, V3, V4 hi cyclic order, the points f(vi), f(v2), /(V3), f(vA) form a conformal 
square. 

The following three points should motivate this definition. 

• Like the definition of isothermic surfaces, this definition of discrete isother- 
mic surfaces is Mobius-invariant. 

• If / : R 2 D D — > R 3 is an immersion, then for e — > 0, 
cr(f(x-e,y-e)J(x + e,y-e),f(x + e,y + e),f(x-e,y + e)) = -l + 0(e 2 ) 
for all x € D if and only if / is an isothermic immersion (see 0). 

• The Christoffel transformation, which also characterizes isothermic sur- 
faces, has a natural discrete analogue (see Propositions ^ an d EJ • The 
condition that all vertex degrees have to be even is used in Proposition [21 

Interior vertices with degree different from 4 play the role of umbilic points. 
At all other vertices, two edge paths — playing the role of curvature lines — 
intersect transversally. It is tempting to visualize a discrete isothermic surface 
as a polyhedral surface with planar quadrilateral faces. However, one should 
keep in mind that those planar faces are not Mobius invariant. On the other 
hand, when we will define discrete minimal surfaces as special discrete isothermic 
surfaces, it will be completely legitimate to view them as polyhedral surfaces 
with planar faces because the class of discrete minimal surfaces is not Mobius 
invariant anyway. 

The Christoffel transformation jS] (see \H>\ for a modern treatment) trans- 
forms an isothermic surface into a dual isothermic surface. It plays a crucial role 
in our considerations. For the reader's convenience, we provide a short proof of 
Proposition n] 

Proposition 1. Let f : R 2 D D — ► R 3 be an isothermic immersion, where D 
is simply connected. Then the formulas 

Ja, ~n/*ir Jy ~ \\fy\\ 2 ' {1 

define (up to a translation) another isothermic immersion f* : R 2 D D —> R 3 
which is called the Christoffel transformed or dual isothermic surface. 
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Proof. First, we need to show that the 1-form df* = f* dx + f* dy is closed and 
thus defines an immersion /*. From equations we have f xy = af x + bf y , 
where a and b are functions of x and y. Taking the derivative of equations (J5J) 
with respect to y and x, respectively, we obtain 

fxy = -|Jj2 (~ a /* + bfy) = - JjTpiaf* - h fv) = fyx- 

Hence, df* is closed. Obviously, ||/*|| = ||/*||, /|_L/*, and f* y G span{/*,/*}. 
Hence, /* is isothermic. □ 

Remarks, (i) In fact, the Christoffel transformation characterizes isothermic 
surfaces: If / is an immersion and equations J5J) do define another surface, then 
/ is isothermic. 

(ii) The Christoffel transformation is not Mobius invariant: The dual of a 
Mobius transformed isothermic surface is not a Mobius transformed dual. 

(Hi) In equations @, there is a minus sign in the equation for /*, but not 
in the equation for /*. This is an arbitrary choice. Also, a different choice of 
conformal curvature line parameters, this means choosing (Ax, Xy) instead of 
{x,y), leads to a scaled dual immersion. Therefore, it makes sense to consider 
the dual isothermic surface as defined only up to translation and (positive or 
negative) scale. 

The Christoffel transformation has a natural analogue in the discrete setting: 
In Proposition |2 we define the dual discrete isothermic surface. The basis for 
the discrete construction is the following Lemma. Its proof is straightforward 
algebra. 

Lemma 1. Suppose a, b, a', b' € C \ {0} with 



a a 



a + b + a' + b' = 0, T77 = -l 
oo' 



id let 



a" h V 



a a b b 

where ~z denotes the complex conjugate of z . Then 

* /* 

a +b +a +b =0, — = -1. 

Proposition 2. Let f : V(T>) — > R 3 be a discrete isothermic surface, where the 
quad-graph T> is simply connected. Then the edges ofT> may be labelled "+" and 
"—" such that each quadrilateral has two opposite edges labelled "+" and the 
other two opposite edges labeled "—" (see Fig. ^p. The dual discrete isothermic 
surface is defined by the formula 

Af* = ± 



A/ir 



where Af denotes the difference of neighboring vertices and the sign is chosen 
according to the edge label. 
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Figure 4: Edge labels of a discrete isothermic surface. 



For a consistent edge labelling to be possible, it is necessary that each vertex 
has an even number of edges. This condition is also sufficient if the the surface 
is simply connected. 

In Definition we define S-quad-graphs. These are specially labeled quad- 
graphs that are used in Definition^of S-isothcrmic surfaces. These S- isothermic 
surfaces form the subclass of discrete isothermic surfaces that we use to define 
discrete minimal surfaces in Section For a discussion why S-isothermic sur- 
faces are the right class to consider, see the Remark at the end of Section 

Definition 3. An S-quad-graph is a quad-graph T> with interior vertices of even 
degree as in Definition and the following additional properties (see Fig. EJ: 

(i) The f-skeleton of T> is bipartite and the vertices are bi-colored "black" and 
"white" . (Then each quadrilateral has two black vertices and two white 
vertices.) 

(ii) Interior black vertices have degree 4. 

(Hi) The white vertices be labeled © and © in such a way that each quadri- 
lateral has one white vertex labeled © and one white vertex labeled ©. 

Definition 4. Let T> be an S-quad-graph, and let V&(X>) be the set of black 
vertices. A discrete S-isothermic surface is a map 

fb : V b (V) E 3 , 

with the following properties: 

(i) If «i, . . . ,v 2n G Vb(T>) are the neighbors of a ©-labeled vertex in cyclic 
order, then /&(i>i), . . . , /&(t>2n) he on a circle in M 3 in the same cyclic order. 
This defines a map from the ©-labeled vertices to the set of circles in K 3 . 

(ii) If vi, . . . ,V2n G Vb(T>) are the neighbors of an ©-labeled vertex, then 
fb(vi), ■ ■ ■ , fb(v2 n ) lie on a sphere in R 3 . This defines a map from the 
©-labeled vertices to the set of spheres in M 3 . 

(Hi) If v c and v s are the ©-labeled and the ©-labeled vertex of a quadrilateral 
of V, then the circle corresponding to v c intersects the sphere correspond- 
ing to v s orthogonally. 

There are two spheres through each black vertex, and the orthogonality 
condition (Hi) of Definition 0] implies that they touch. Likewise, the two circles 
at a black vertex touch, i.e. they have a common tangent at the single point of 
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Figure 5: Left: Schramm's circle patterns as discrete conformal maps. Right: The 
combinatorics of S-quad-graphs. 

intersection. Discrete S-isothermic surfaces are therefore composed of touching 
spheres and touching circles with spheres and circles intersecting orthogonally. 
Interior white vertices of degree unequal to 4 are analogous to umbilic points 
of smooth isothermic surfaces. Generically, the orthogonality condition (Hi) 
follows from the seemingly weaker condition that the two circles through a 
black vertex touch: 

Lemma 2. (Touching Coins Lemma) Whenever four circles in 3- space touch 
cyclically but do not lie on a common sphere, they intersect the sphere which 
passes through the points of contact orthogonally. 

From any discrete S-isothermic surface, one obtains a discrete isothermic 
surface (as in Definitional by adding the centers of the spheres and circles: 

Definition 5. Let /& : Vb(T>) — > K 3 be a discrete S-isothermic surface. The 
central extension of fb is the discrete isothermic surface / : V — > K 3 defined by 

f(v) = f b (v) ifveV b , 

and otherwise 

f(v) = the center of the circle or sphere corresponding to v. 

The central extension of a discrete S-isothermic surface is indeed a discrete 
isothermic surface: The quadrilaterals corresponding to the faces of the quad- 
graph are planar right-angled kites (see Fig. (right)) and therefore conformal 
squares. 

The following statement is easy to see 0]. It says that the duality transfor- 
mation preserves the class of discrete S-isothermic surfaces. 

Proposition 3. The Christoffel dual of a central extension of a discrete S- 
isothermic surface is a central extension of a discrete S-isothermic surface. 

The construction of the central extension does depend on the choice of a 
point at infinity, because the centers of circles and spheres are not Mobius in- 
variant. Strictly speaking, a discrete S-isothermic surface has a 3-paramcter 
family of central extensions. However, we will assume that one infinite point 
is chosen once and for all and we will not distinguish between S-isothermic 
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Q touching spheres 

i / orthogonal circles 

| | planar faces 

P~] orthogonal kite 



Figure 6: Geometry of a discrete S-isothcrmic surface without "umbilics" . 

surfaces and their central extension. Then it also makes sense to consider the 
S-isothermic surfaces as polyhedral surfaces. Note that all planar kites around a 
©-labeled vertex lie in the same plane: the plane that contains the correspond- 
ing circle. We will therefore consider an S-isothermic surface as a polyhedral 
surface whose faces correspond to ©-labeled vertices of the quad-graph, whose 
vertices correspond to ©-labeled vertices of the quad-graph, and whose edges 
correspond to the black vertices of the quad-graph. The elements of a discrete 
S-isothcrmic surface are shown schematically in Fig. El Hence: 

A discrete S-isothermic surface is a polyhedral surface such that the faces 
have inscribed circles and the inscribed circles of neighboring faces touch their 
common edge in the same point. 

In view of the Touching Coins Lemma (Lemma |5J), this could almost be an 
alternative definition. 

The following lemma, which follows directly from Lemma E] describes the 
dual discrete S-isothcrmic surface in terms of the corresponding polyhedral dis- 
crete S-isothermic surface. 

Lemma 3. Let P be a planar polygon with an even number of cyclically ordered 
edges given by the vectors 1%, . . ., ?2n G K 2 , h + - ■ - + hn = 0. Suppose the polygon 
has an inscribed circle with radius R. Let rj be the distances from the vertices 
of P to the nearest touching point on the circle: \\lj\\ = rj + ry+i. Then the 
vectors I*, . . . , l% n given by 

I* = {-ly-t—ii 

form a planar polygon with an inscribed circle with radius 1/R. 

It follows that the radii of corresponding spheres and circles of a discrete 
S-isothcrmic surface and its dual are reciprocal. 

3 Koebe polyhedra 

In this section we construct special discrete S-isothermic surfaces, which we call 
the Koebe polyhedra, coming from circle packings (and more general orthogonal 
circle patterns) in S 2 . 
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Figure 7: Left: A circle packing corresponding to a triangulation. Middle: 
The orthogonal circles. Right: A circle packing corresponding to a cellular 
decomposition with orthogonal circles. 

A circle packing in S 2 is a configuration of disjoint discs which may touch 
but not intersect. Associating vertices to the discs and connecting the vertices 
of touching discs by edges one obtains a combinatorial representation of a circle 
packing, see Fig. (left) . 

In 1936, Koebe published the following remarkable statement about circle 
packings in the sphere |22j . 

Theorem 1. (Koebe) For every triangulation of the sphere there is a packing 
of circles in the sphere such that circles correspond to vertices, and two circles 
touch if and only if the corresponding vertices are adjacent. This circle packing 
is unique up to Mobius transformations of the sphere. 

Observe that for a triangulation one au- 
tomatically obtains not one but two orthog- 
onally intersecting circle packings as shown 
in Fig. [3 (middle). Indeed, the circles pass- 
ing through the points of contact of three 
mutually touching circles intersect these or- 
thogonally. This observation leads to the 
following generalization of Koebe's theorem Figure 8" The Koebe polyhedron 
to cellular decompositions of the sphere with as a discrete S-isothermic surface, 
faces which are not necessarily triangular, 
see Fig. (right). 

Theorem 2. For every polytopal 1 cellular decomposition of the sphere, there 
exists a pattern of circles in the sphere with the following properties. There is 
a circle corresponding to each face and to each vertex. The vertex circles form 
a packing with two circles touching if and only if the corresponding vertices are 
adjacent. Likewise, the face circles form a packing with circles touching if and 
only if the corresponding faces are adjacent. For each edge, there is a pair of 
touching vertex circles and a pair of touching face circles. These pairs touch in 
the same point, intersecting each other orthogonally. 

This circle pattern is unique up to Mobius transformations. 

The first published statement and proof of this theorem seems to be con- 
tained in 0. For generalizations, see [211, an( i P3 > the latter also for a 
variational proof (see also Section |H1 of this article) . 

Now, mark the centers of the circles with white dots and mark the intersec- 
tion points, where two touching pairs of circles intersect each other orthogonally, 
with black dots. Draw edges from the center of each circle to the intersection 

1 We call a cellular decomposition of a surface polytopal, if the closed cells are closed discs, 
and two closed cells intersect in one closed cell if at all. 
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points on its periphery. You obtain a quad-graph with bicolored vertices. Since, 
furthermore, the black vertices have degree four, the white vertices may be 
labeled ® and © to make the quad-graph an S-quad-graph. 

Now let us construct the spheres intersecting S 2 orthogonally along the cir- 
cles marked by ®. Connecting the centers of touching spheres, one obtains a 
Koe.be polyhedron: a convex polyhedron with all edges tangent to the sphere S 2 . 
Moreover, the circles marked with © are inscribed into the faces of the polyhe- 
dron, sec Fig. IS] Thus we have a polyhedral discrete S-isothermic surface. The 
discrete S-isothcrmic surface is given by the spheres @ and the circles ©. 

Thus, Theorem implies the following theorem. 

Theorem 3. Every polytopal cell decomposition of the sphere can be realized by 
a polyhedron with edges tangent to the sphere. This realization is unique up to 
projective transformations which fix the sphere. 

There is a simultaneous realization of the dual polyhedron, such that corre- 
sponding edges of the dual and the original polyhedron touch the sphere in the 
same points and intersect orthogonally. 

The last statement of the theorem follows from the construction if one in- 
terchanges the © and © labels. 



4 Discrete minimal surfaces 

The following theorem about continuous minimal surfaces is due to Christof- 
fel [Sj. For a modern treatment, see |15|. This theorem is the basis for our 
definition of discrete minimal surfaces. We provide a short proof for the reader's 
convenience. 

Theorem 4 (Christoffel). Minimal surfaces are isothermic. An isothermic 
immersion is a minimal surface, if and and only if the dual immersion is con- 
tained in a sphere. In that case the dual immersion is in fact the Gauss map of 
the minimal surface, up to scale and translation. 

Proof. Let / be an isothermic immersion with normal map N. Then 

(N x , f x ) = X 2 ki and (N y , f y ) = X 2 k 2 , 

where k\ and ki are the principal curvature functions of / and A = 1 1 f x \ \ = \\f y \\. 
By equations (0), the dual isothermic immersion /* has normal N* = —N, and 



(Kj;) = {-N y , 



fy 



Its principal curvature functions are therefore 

k 1 = — — and fc 2 = 

Hence / is minimal (this means k\ = —kz) if and only if /* is contained in 
a sphere (k* = k%). In that case, /* is the Gauss map of / (up to scale and 
translation), because the tangent planes of / and /* at corresponding points 
are parallel. □ 
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The idea is to define discrete minimal surfaces as S-isothcrmic surfaces which 
are dual to Koebe polyhcdra; the latter being a discrete analogue of conformal 
paramctcrizations of the sphere. By thcorcmElbclow, this leads to the following 
definition. 

Definition 6. A discrete minimal surface is an S-isothermic discrete surface 
F : Q — ► R 3 which satisfies any one of the equivalent conditions (i)-(iii) below. 

Suppose x G Q is a white vertex of the quad-graph Q such that F(x) is the 
center of a sphere. Let y\ . . . ym be the vertices neighboring x in Q. (Generically, 
ri = 2.) Then F(yj) are the points of contact with the neighboring spheres and 
simultaneously points of intersection with the orthogonal circles. Let F(yj) = 
F(x) + bj. (Sec figure^]) Then the following equivalent conditions hold: 

(i) The points F(x) + (— l) J 6j lie on a circle. 

(ii) There is an TV G 1R 3 such that (— 1) J (bj, N) is the same for j = 1, . . . , In. 

(iii) There is plane through F(x) such that the points {F(yj) | j even} and the 
points {F(yj) \ j odd} lie in planes which arc parallel to it at the same 
distance on opposite sides. 

Remark. The definition implies that a discrete minimal surface is a polyhedral 
surface with the property that every interior vertex lies in the convex hull of its 
neighbors. This is the maximum principle for discrete minimal surfaces. 

Examples. Fig. ^ (left) shows a discrete minimal Enneper surface. Only the 
circles arc shown. A variant of the discrete minimal Enneper surface is shown 
in Fig. 1161 Here, only the touching spheres are shown. Fig. ^ (right) shows a 
discrete minimal catenoid. Both spheres and circles are shown. Fig. |21 shows a 
discrete minimal Schwarz P-surface and a discrete minimal Scherk tower. 

These examples are discussed in detail in section ITTTl 

Theorem 5. An S-isothermic discrete surface is a discrete minimal surface, if 
and only if the dual S-isothermic surface corresponds to a Koebe polyhedron. 

Proof. That the S-isothcrmic dual of a Koebe polyhedron is a discrete minimal 
surface is fairly obvious. On the other hand, let F : Q — > R 3 be a discrete 
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minimal surface and let x € Q and yi ■ ■ ■ ym G Q be as in Definition Let 
F : Q — >JR 3 be the dual S-isothermic surface. We need to show that all 
circles of F lie in one and the same sphere S and that all the spheres of F 
intersect S orthogonally. It follows immediately from DefmitionElthat the points 
F(yi) . . . F(y 2n ) he on a circle c x in a sphere S x around F(x). Let S be the 
sphere which intersects S x orthogonally in c x . The orthogonal circles through 
F(yi) . . . F(y 2 n) also lie in S. Hence, all spheres of F intersect S orthogonally 
and all circles of F lie in S. □ 

Remark. Why do we use S-isothcrmic surfaces to define discrete minimal sur- 
faces? Alternatively, one could define discrete minimal surfaces as the surfaces 
obtained by dualizing discrete (cross-ratio —1) isothermic surfaces with all quad- 
graph vertices in a sphere. Indeed, this definition was proposed in However, 
it turns out that the class of discrete isothermic surfaces is too general to lead 
to a satisfactory theory of discrete minimal surfaces. 

Every way to define the concept of a discrete isothermic immersion imposes 
an accompanied definition of discrete conformal maps. Since a conformal map 
R 2 D D — > R 2 is just an isothermic immersion into the plane, discrete confor- 
mal maps should be defined as discrete isothermic surfaces that lie in a plane. 
Definition |21 for isothermic surfaces implies the following definition for discrete 
conformal maps: A discrete conformal map is a map from a domain of 7Z? to the 
plane such that all elementary quads have cross-ratio —1. The so defined dis- 
crete conformal maps are too flexible. In particular, one can fix one sublattice 
containing every other point and vary the other one, see 0]. 

Definition 0] for S-isothermic surfaces, on the other hand, leads to discrete 
conformal maps that are Schramm's "circle patterns with the combinatorics of 
the square grid" |26j . This definition of discrete conformal maps has many ad- 
vantages: First, there is Schramm's convergence result (ibid). Secondly, orthog- 
onal circle patterns have the right degree of rigidity. For example, by Theorem|2 
two circle patterns that correspond to the same quad-graph decomposition of 
the sphere differ by a Mobius transformation. One could say: The only discrete 
conformal maps from the sphere to itself are the Mobius transformations. Fi- 
nally, a conformal map / : R 2 D D — > R 2 is characterized by the conditions 

\M = \fy\, fx±fy (3) 

To define discrete conformal maps / : 1? D D — > C, it is natural to impose 
these two conditions on two different sub-lattices (white and black) of Z 2 , i.e. to 
require that the edges meeting at a white vertex have equal length and the edges 
at a black vertex meet orthogonally. Then the elementary quadrilaterals are 
orthogonal kites, and discrete conformal maps are therefore precisely Schramm's 
orthogonal circle patterns. 

5 A Weierstrass-type representation 

In the classical theory of minimal surfaces, the Weierstrass representation allows 
the construction of an arbitrary minimal surface from holomorphic data on the 
underlying Ricmann surface. We will now derive a formula for discrete minimal 
surfaces that resembles the Weierstrass representation formula. An orthogonal 
circle pattern in the plane plays the role of the holomorphic data. The discrete 
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Weierstrass representation describes the S-isothermic minimal surface that is 
obtained by projecting the pattern stereographically to the sphere and dualizing 
the corresponding Koebe polyhedron. 



Theorem 6 (Weierstrass representation). Let Q be an S- quad- graph, and 
let c : Q — > (D be an orthogonal circle pattern in the plane: For white vertices 
x G Q, c(x) is the center of the corresponding circle, and for black vertices 
y G Q, c(y) is the corresponding intersection point. The S-isothermic minimal 
surface 

F : {x G Q | x is labelled ® } — > R 3 , 
F(x) = the center of the sphere corresponding to x 

that corresponds to this circle pattern is given by the following formula. Let 
xi,%2 G Q be two vertices, both labelled (s), that correspond to touching circles 
of the pattern, and let y G Q be the black vertex between x\ and x 2 , which 
corresponds to the point of contact. The centers F{x\) and F(x,2) of the corre- 
sponding touching spheres of the S-isothermic minimal surface F satisfy 



F{x 2 )-F( Xl ) 



±Rc 



R(x 2 ) + R{x\) c(x 2 ) - c{x\) 



1 + H 2 \c(x 2 ) - c(xi)\ 
where p = c(y) and the radii R(xj) of the spheres are 
R( Xj ) 




(4) 



l + | c(a;j .)|2-|cfe)-jf 



2|c(^)-p| 



(5) 



The sign on the right hand side of equation depends on whether the two 
edges of the quad-graph connecting x\ with y and y with x 2 are labelled '+' or 
'-' (see Figs. g]cmd[J (right)/ 

Proof. Let s : <C — ► S 2 C R 3 be the stereographic projection 



s(p) 



2 Rep 
2 Imp 

>| 2 -1, 



Its differential is 



ds p (v) 




and 



|ds P (u)| 



2\v\ 



l + |pP' 



where || • || denotes the Euclidean norm. 

The edge between F(xi) and F(x 2 ) of F has length Ri + R 2 (this is obvi- 
ous) and is parallel to ds v (c{x 2 ) — c(xi)). Indeed, this edge is parallel to the 
corresponding edge of the Koebe polyhedron, which, in turn, is tangential to 
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the orthogonal circles in the unit sphere, touching in c{p). The pre- images of 
these circles in the plane touch in p, and the contact direction is c(x2) — c(x\). 
Hence, equation 0} follows from 

/ * , ^ / , v , a ds v (c(x2) — c(Xi)) 

F(x 2 ) - F( Xl ) = ±{R(x 2 ) + R( Xl )) p ) . 

\\ds p {c{X2) - c{xi))\\ 

To show equation |J5J, note that the stereographic projection s is the restric- 
tion of the reflection on the sphere around the north pole N of S 2 with radius 
y/2, restricted to the equatorial plane <D. See Fig. QUI We denote this reflection 
also by s. Consider the sphere with center c(xj) and radius r = \c(xj) — p\, 
which intersects the equatorial plane orthogonally in the circle of the planar 
pattern corresponding to Xj . This sphere intersects the ray from the north pole 
N through c(xj) orthogonally at the distances d±r from N, where d, the dis- 
tance between N and c(xj), satisfies d 2 = 1 + |c(xj)| 2 . This sphere is mapped 
by s to a sphere, which belongs to the Koebe polyhedron and has radius 1/Rj. 
It intersects the ray orthogonally at the distances 2/(d±r). Hence, its radius is 



2r 



d 2 -7 

Equation JSJ) follows. □ 



6 The associated family 

Every continuous minimal surfaces comes with an associated family of isometric 
minimal surfaces with the same Gauss map. Catenoid and helicoid are members 
of the same associated family of minimal surfaces. The concept of an associ- 
ated family carries over to discrete minimal surfaces. In the smooth case, the 
members of the associated family remain conformally, but not isothermically, 
parameterized. Similarly, in the discrete case, one obtains discrete surfaces 
which are not S-isothermic but should be considered as discrete conformally 
parameterized minimal surfaces. 

The associated family of an S-isothermic minimal surface consists of the 
one-parameter family of discrete surfaces that are obtained by the following 
construction. Before dualizing the Koebe-polyhedron (which would yield the 
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Figure 11: Proof of Lemma|H The vector Vj is obtained by rotating Vj in the 
tangent plane to the sphere at Cj. 

S-isothcrmic minimal surface), rotate each edge by an equal angle in the plane 
which is tangent to the unit sphere in the point where the edge touches the unit 
sphere. 

This construction leads to well defined surfaces because of the following 
lemma, which is an extension of Lemma Sec Fig. 1111 

Lemma 4. Let P be a planar polygon with an even number of cyclically ordered 
edges given by the vectors l\, . . . , li n € 1R 3 ,, + . . . + ?2?i = 0. Suppose the polygon 
has an inscribed circle c with radius R, which lies in a sphere S. Let rj be the 
distances from the vertices of P to the nearest touching point on the circle: 
\\lj\\ = rj + r J+ i. Rotate each vector lj by an equal angle tp in the plane which is 
tangent to S in the point Cj where the edge touches the S to obtain the vectors 
l\ , . . . , 4n ■ Then the vectors i[ , . . . , 4« given by 

= t_ l y_±_ l {v) 

satisfy l^* + . . . + l^* = 0; that is, they form a (non-planar) closed polygon. 

Proof. For j = 1, . . . , 2n let (vj, Wj, n) be the orthonormal basis of R 3 which is 
formed by Vj = lj/\\lj\\, the unit normal n to the plane of the polygon P, and 

Wj = nxvj. (6) 

Let Vj be the vector in the tangent plane to the sphere S at Cj that makes an 
angle (p with Vj. Then 

Vj = cos ip Vj + sin ip cos 8 Wj + sin ip sin 8 n, 

where 8 is the angle between the tangent plane and the plane of the polygon. 
This angle is the same for all edges. Since 



(v) 
j ' 
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we have to show that 



^^(— l) 3 ( 1 j (cosies Vj + sin ^ cos # Wj + sin ip sin 6 n) = 0. 

By Lemma 

2n , 



i=i 



Due to ©, 



as well. Finally, 



2 ™ / 1 



j'=i 



//'; 



2n , 



= 0, 



because it is a telescopic sum. □ 

The following two theorems are easy to prove. First, the Weierstrass-type 
formula of Theorem may be extended to the associate family. 

Theorem 7. Using the notation of Theorem^ the discrete surfaces F v of the 
associated family satisfy 




R(x 2 ) + R(xi) c(x 2 ) - c{xi) 




1 + H 2 \c{x 2 ) - c(xi)| 

Fig. 1 121 shows the associated family of the S-isothermic catenoid. The essen- 
tial properties of the associated family of a continuous minimal surface — that 
the surfaces are isometric and have the same Gauss map — carries over to the 
discrete setting in the following form. 

Theorem 8. The surfaces F v of the associated family of an S-isothermic min- 
imal surface Fq consist, like Fq, of touching spheres. The radii of the spheres 
do not depend on tp. 

In the generic case, when the quad-graph has H? -combinatorics, there are 
also circles through the points of contact, like it is the case with Fq . The normals 
of the circles do not depend on tp. 

This theorem follows directly from the geometric construction of the associ- 
ated family (Lemma QJ. 



7 Convergence 

Schramm has proved the convergence of circle patterns with the combinatorics 
of the square grid to meromorphic functions [23 • Together with the Weierstrass- 
type representation formula for S-isothermic minimal surfaces, this implies the 
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Figure 12: The associated family of the S-isothermic catenoid. The Gauss map 
is preserved 

following approximation theorem for discrete minimal surfaces. Fig. 1131 illus- 
trates the convergence of S-isothermic Enneper surfaces to the continuous En- 
neper surface. 

Theorem 9. Let D C (D be a simply connected bounded domain with smooth 
boundary, and let W C <C be an open set that contains the closure of D. Suppose 
that F : W — ► 1R 3 is a minimal immersion without umbilic points in conformal 
curvature line coordinates. There exists a sequence of S-isothermic minimal 
surfaces F n : Q n — > M 3 such that the following holds. Each Q n is a simply 
connected S-quad-graph in D which is a subset of the lattice \T1? ■ If for x G D, 
we define F n (x) to be the value of F n at a point of Q n closest to x, then F n 
converges to F uniformly with error 0(— ) on compacts in D. In fact, the whole 
associated families F UtV converge to the associated family F v of F uniformly 
(also in (p) and with error 0(— ) on compacts in D. 

Proof. Assuming that F is appropriately scaled, 




where g : W — ► (D is a locally injective meromorphic function. By Schramm's 
results (Theorem 9.1 of and the remark on p. 387), there exists a sequence of 
orthogonal circle patterns c„ : Q n — » (D approximating g and g' uniformly and 
with error 0(— ) on compacts in D. Define F n by the Weierstrass formula (J3J 
with data c — c n . Using the notation of Theorem one finds 

! j / 1-S(V) 2 \ ! , ! . 
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Figure 13: A sequence of S-isothermic minimal Enneper surfaces in different 
discretizations. 

uniformly on compacts in D. Setting F n = -^F n , the convergence claim follows. 
The same reasoning applies to the whole associated family of F. □ 



8 Orthogonal circle patterns in the sphere 

In the simplest cases, like the discrete Enneper surface and the discrete catenoid 
(Fig. ^) , the construction of the corresponding circle patterns in the sphere can 
be achieved by elementary methods, see Section ITUl In general, the problem is 
not elementary. Developing methods introduced by Colin de Verdiere jH], the 
first and third author have given a constructive proof of the generalized Koebe 
theorem, which uses a variational principle pj- It also provides a method for 
the numerical construction of circle patterns (see also [27]). An alternative 
algorithm was implemented in Stephenson's program circlepack 1 2 . It is 
based on methods developed by Thurston |29|. The first step in both methods 
is to transfer the problem from the sphere to the plane by a stereographic 
projection. Then the radii of the circles are calculated. If the radii are known, 
it is easy to reconstruct the circle pattern. The radii are determined by a set 
of nonlinear equations, and the two methods differ in the way in which these 
equations are solved. Thurston-type methods work by iteratively adjusting the 
radius of each circle so that the neighboring circles fit around. The above 
mentioned variational method is based on the observation that the equations 
for the radii are the equations for a critical point of a convex function of the 
radii. The variational method involves minimizing this function to solve the 
equations. 

Both of these methods may be used to construct the circle patterns for 
the discrete Schwarz P-surfacc and for the discrete Scherk tower, see Figs. 
and 1151 One may also take advantage of the symmetries of the circle patterns 
and construct only a piece of it (after stereographic projection) as shown in 
Fig. To this end, one solves the Euclidean circle pattern problem with 
Neumann boundary conditions: For boundary circles, the nominal angle to be 
covered by the neighboring circles is prescribed. 

However, we have actually constructed the circle patterns for the discrete 
Schwarz P-surface and the discrete Scherk tower using an new method suggested 
in [2B1 • It is a variational method that works directly on the sphere. No stereo- 
graphic projection is necessary; the spherical radii of the circles are calculated 
directly. The variational principle for spherical circle patterns is completely 
analogous to the variational principles for Euclidean and hyperbolic patterns 
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Figure 14: A piece of the circle pattern for a Schwarz P-surface after stereo- 
graphic projection to the plane. 

presented in pj. We briefly describe our variational method for circle patterns 
on the sphere. For a more detailed exposition, the reader is referred to |28j . 
Here, we will only treat the case of orthogonally intersecting circles. 

The spherical radius r of a non-degenerate circle in the unit sphere satisfies 
< r < ir. Instead of the radii r of the circles, we use the variables 

p = logtan(r/2). (8) 

For each circle j, we need to find a pj such that the corresponding radii solve 
the circle pattern problem. 

Proposition 4. The radii Tj are the correct radii for the circle pattern if and 
only if the corresponding pj satisfy the following equations, one for each circle: 
The equation for circle j is 

2^(arctane pfc_ft + arctane Pfe+ft ) = (9) 

neighbors k 

where the sum is taken over all neighboring circles k. For each circle j , <5j is the 
nominal angle covered by the neighboring circles. It is normally 2tt for interior 
circles, but it differs for circles on the boundary or for circles where the pattern 
branches. 

The equations © are the equations for a critical point of the functional 
S(p) = (lmLi 2 (ie pfc ~ ft ) +ImLi 2 (ie«- pfc ) 

3 
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Here, the first sum is taken over all pairs (j, k) of neighboring circles, the second 
sum is taken over all circles j. The diloqarithm function Li2(z) is defined by 

Li 2 (z) = -pog(l-C«C- 

The second derivative of S(p) is the quadratic form 

D 2 S = V ( — (dp k - d Pj f — (dp k + d Pi ) 2 ) , (10) 

where the sum is taken over pairs of neighboring circles. 

We provide a proof of Proposition 0] in the Appendix. Unlike the analogous 
functionals for cuclidcan and hyperbolic circle patterns, the functional S( P ) is 
unfortunately not convex : The second derivative is negative for the tangent 
vector v = ^2jd/d P j. The index is therefore at least 1. Thus, one cannot 
simply minimize S to get to a critical point. However, the following method 
seems to work. Define a reduced functional S(p) by maximizing in the direction 
v. 

S(p) = umx.S(p + tv). (11) 

Obviously, S( P ) is invariant under translations in the direction v. Now the idea 
is to minimize S( P ) restricted to £^ . P j = 0. This method has proved to be 
amazingly powerful. In particular, it can be used to produce branched circle 
patterns in the sphere. It would be very interesting and important to give a 
theoretical explanation of this phenomenon. 

9 Constructing discrete minimal surfaces 

Given a specific continuous minimal surface, how does one construct an analo- 
gous discrete minimal surface? In this section we outline the general method for 
doing this. By way of an example, Figure IT51 illustrates the construction of the 
Schwarz P-surfacc. Details on how to construct the concrete examples in this 
paper are explained in Section ITU1 The difficult part is finding the right circle 
pattern (paragraphs 1 and 2 below). The remaining steps, building the Koebe 
polyhedron and dualizing it (Paragraphs 3 and 4) are rather mechanical. 

1. Investigate the Gauss image of the curvature lines. The Gauss map 
of the continuous minimal surface maps its curvature lines to the sphere. Obtain 
a qualitative picture of this image of the curvature lines under the Gauss map. 
This yields a quad-graph immersed in the sphere. Here one has to choose how 
many curvature lines one wants to use. This corresponds to a choice of different 
levels of refinement of the discrete surface. Also, a choice is made as to which 
vertices will be black and which will be white. This choice is usually determined 
by the nature of the exceptional vertices corresponding to umbilics and ends. 
(Umbilics have to be black vertices.) Only the combinatorics of this quad graph 
matter. (Figure top left.) Generically, the (interior) vertices have degree 4. 
Exceptional vertices correspond to ends and umbilic points of the continuous 
minimal surface. In the Figurc ll5l the corners of the cube are exceptional. They 
correspond to the umbilic points of the Schwarz P-surfacc. The exceptional 
vertices may need to be treated specially. For details see section ITU1 
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Gauss image of the curvature lines 



circle pattern 




Koebe polyhedron discrete minimal surface 



Figure 15: Construction of the discrete Schwarz P-surface. 

2. Construct the circle pattern. From the quad graph, construct the corre- 
sponding circle pattern. White vertices will correspond to circles, black vertices 
to intersection points. Usually, the generalized Koebe theorem is evoked to 
assert existence and Mobius uniqueness of the pattern. The problem of practi- 
cally calculating the circle pattern was discussed in Section 00 Use symmetries 
of the surface or special points where you know the direction of the normal to 
eliminate the Mobius ambiguity of the circle pattern. 

3. Construct the Koebe polyhedron. From the circle pattern, construct 
the Koebe polyhedron. Here, a choice is made as to which circles will become 
spheres and which will become circles. The two choices lead to different discrete 
surfaces close to each other. Both are discrete analogues of the continuous 
minimal surface. 

4. Discrete minimal surface. Dualize the Koebe polyhedron to obtain a 
minimal surface. 

If the function g(z) in the Weierstrass representation (JZJ) of the continuous 
minimal surface is simple enough, it may be that one can construct an orthogonal 
circle pattern that is analogous to (or even approximates) this holomorphic 
function explicitly by some other means. For example, this is the case for the 
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Enneper surfaces and the catenoid (see Section[5J). In this case one does not use 
Koebe's theorem to construct the circle pattern. 



10 Examples 

We now apply the method outlined in the previous section to construct concrete 
examples of discrete minimal surfaces. In the case of Enneper's surface, the 
orthogonal circle pattern is trivial. The circle patterns for the higher order 
Enneper surfaces and for the catenoid are known circle pattern analogues of the 
functions z a and e z . To construct the circle patterns for the Schwarz P-surface 
and the Scherk tower, we use Koebe's theorem. 



10.1 Enneper's surface 



The Weierstrass representation of Enneper's surface in conformal curvature line 
coordinates is equation J7J) with g(z) = z. The domain is (D, and there are no 
umbilic points. In the domain, the curvature lines arc the parallels to the real 
and imaginary axes. The Gauss map embeds the domain into the sphere. 

The quad graph that captures this qual- 
itative behavior of the curvature lines is the 
regular square grid decomposition of the plane. 
There are also obvious candidates for the 
circle patterns to use: Take an infinite reg- 
ular square grid pattern in the plane. It 
consists of circles with equal radius r and 
centers on a square grid with spacing \/2r. 
It was shown by He ^31 that these pat- 
terns are the only embedded and locally 
finite orthogonal circle patterns with this 
quad graph. Project it stereographically to 
the sphere and build the Kocbc polyhedron. 
Dualize to obtain a discrete version of En- 
neper's surface. See Figs.Q] (left) and El 




Figure 16: S-isothermic higher 
order Enneper surface. Only the 
spheres are shown. 



10.2 The higher order Enneper surfaces 

As the next example, consider the higher order Enneper surfaces Their 
Weierstrass representation has g(z) = z a . One may think of them as Enneper 
surfaces with an umbilic point in the center. 

An orthogonal circle pattern analogue of the maps z a was introduced in 0] . 
Sectors of these circle patterns were proven to be embedded |2j, CO- Stereo- 
graphic projection to the sphere followed by dualization leads to S-isothermic 
analogues of the higher order Enneper surfaces. An S-isothermic higher order 
Enneper surface with a simple umbilic point (a = 4/3) is shown in Fig. 1161 



10.3 The catenoid 

The next most simple example is a discrete version of the catenoid. Here, 
g(z) = e z . The corresponding circle pattern in the plane is the S-Exp pattern 0], 



24 



a discretization of the exponential map. The underlying quad-graph is U? , with 
circles corresponding to points (m,n) with to + n = mod 2. The centers 
c(to, n) and the radii r(n, to) of the circles are 

c(n,m) = e an+ip "\ r(n,m) = sm(p)\c(n, m)|, 

where 

p = ir/N, a = arctanh(i|l - e 2ip |). 

(It is not true that c(m, n) is an intersection point if m + n = 1 mod 2.) 

The corresponding S-isothermic minimal surface is shown in Fig.^ (right). 
The associated family of the discrete catenoid (see SectionlHJ is shown in Fie . 1121 

Other discrete versions of the catenoid have been put forward. A discrete 
isothermic catenoid is constructed in This construction can be generalized 
in such a way that one obtains the discrete S-isothermic catenoid described 
above. This works only because the surface is so particularly simple. Then, 
there is also the discrete catenoid constructed in |23| . It is an area minimizing 
polyhedral surface. This catenoid is not related to the S-isothermic catenoid. 

10.4 The Schwarz P-surface 

The Schwarz P-surface is a triply periodic minimal surface. It is the symmetric 
case in a 2-parameter family of minimal surfaces with 3 different hole sizes 
(only the ratios of the hole sizes matter), see [IT]. The domain of the Schwarz 
P-surface, where the translation periods are factored out, is a Riemann surface 
of genus 3. The Gauss map is a double cover of the sphere with 8 branch 
points. The image of the curvature line pattern under the Gauss map is shown 
schematically in Fig.^lffop left), thin lines. It is a refined cube. More generally, 
one may consider three different numbers to, n, and k of slices in the three 
directions. The 8 corners of the cube correspond to the branch points of the 
Gauss map. Hence, not 3 but 6 edges are incident with each corner vertex. The 
corner vertices are assigned the label ©. We assume that the numbers to, n, 
and k are even, so that the vertices of the quad graph may be labelled '©', '©', 
and '•' consistently (see Section |5J). 

To invoke Koebe's theorem (in the form of Theorem forget momentarily 
that we are dealing with a double cover of the sphere. Koebe's theorem implies 
the existence and Mobius-uniqueness of a circle pattern as shown in Fig. 1151 
(top right). (Only one eighth of the complete spherical pattern is shown.) The 
Mobius ambiguity is eliminated by imposing octahedral metric symmetry. 

Now lift the circle pattern to the branched cover, construct the Koebe poly- 
hedron and dualize it to obtain the Schwarz P-surface; see Fig. 1151 (bottom row). 
A fundamental piece of the surface is shown in Fig.0 (left). 

We summarize these results in a theorem. 

Theorem 10. Given three even positive integers m, n, k, there exists a corre- 
sponding unique (unsymmetric) S-isothermic Schwarz P-surface. 

Surfaces with the same ratios to : n : k are different discretizations of the 
same continuous Schwarz P-surface. The cases with to = n = k correspond to 
the symmetric Schwarz P-surface. 
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Figure 17: The combinatorics of the Scherk tower. 
10.5 The Scherk tower 

Finally, consider Scherk's saddle tower, a simply periodic minimal surface, which 
is asymptotic to two intersecting planes. There is a 1-parameter family, the 
parameter corresponding to the angle between the asymptotic planes, see |1 lj . 
An S-isothcrmic minimal Scherk tower is shown in Fig. El (right). 

When mapped to the sphere by the Gauss map, the curvature lines of the 
Scherk tower form a pattern with four special points, which correspond to the 
four half-planar ends. A loop around a special point corresponds to a period of 
the surface. In a neighborhood of each special point, the pattern of curvature 
lines behaves like the image of the standard coordinate net under the map 
z i— > z 2 around z = 0. In the discrete setting, the special points are modeled by 
pairs of 3-valent vertices; see Fig. El fleft)- This is motivated by the discrete 
version of z 2 in 2 . The quad graph we use to construct the Scherk tower 
looks like the quad graph for an unsymmetric Schwarz P-surface with one of the 
discrete parameters equal to 2. The ratio m : n corresponds to the parameter 
of the smooth case. Again, by Koebe's theorem, there exists a corresponding 
circle pattern, which is made unique by imposing metric octahedral symmetry. 
But now we interpret the special vertices differently. Here, they are not branch 
points. The right hand side of Fig.lTTIshows how they are to be treated: Split the 
vertex (and edges) between each pair of 3-valent vertices in two. Then introduce 
new 2-valent vertices between the doubled vertices. Thus, instead of pairs of 
3-valent vertices we now have 2-valent vertices. The newly inserted edges have 
length 0. Thus, stretching the concept a little bit, one obtains infinite edges 
after dualization. This is in line with the fact that the special points correspond 
to half-planar ends. 

Fig. El (right) shows an S-isothermic Scherk tower. 

Theorem 11. Given two even positive integers m and n there exists a corre- 
sponding unique S-isothermic Scherk tower. 

The cases with m = n correspond to the most symmetric Scherk tower, the 
asymptotic planes of which intersect orthogonally. 

Appendix. Proof of Proposition |4] 

Figure E| shows a "flower" of an orthogonal circle pattern: a central circle and 
its orthogonally intersecting neighbors. For simplicity, it shows a circle pattern 
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Figure 18: A flower of an orthogonal circle pattern. 



in the euclidean plane. We are, however, concerned with circle patterns in the 
sphere, where the centers are spherical centers, the radii are spherical radii and 
so forth. The radii of the circles are correct if and only if for each circle the 
neighboring circles "fit around". This means that for each circle j, 

neighbors k 

where tpjk is half the angle covered by circle k as seen from the center of circle j, 
and where normally = 2tt except if j is a boundary circle or a circle where 
branching occurs. (In those cases, $>j has some other given value.) Equations 
follow from the following spherical trigonometry lemma: 

Lemma 5. In a right-angled spherical triangle, let r\ and r^ be the sides en- 
closing the right angle, and let ip be the angle opposite side r^. Then 

tp = arctane p2 ~ Pl + arctan e P2+pl , (12) 

where r and p are related by equation JSJ . 

Proof. Napier's rule says 2 

tanr2 

tan tp = — . 

sm r\ 

Equation 1|12|) follows this and the trigonometric identity 
' tan r2 1 



f tj0il1 To \ 

arctan ( ) = arctane P2_Pl + arctan e P2+Pl . (13) 

V sm n / 

(To derive equation (|13() . start by applying the identity 

a + b 

arctan a + arctan o = arctan ■ 



I -ah 

to its right hand side.) □ 

2 In several editions of Bronshtein and Semendyayev's Handbook of Mathematics there is a 
misprint in the corresponding equation. 
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Now let 

f(x) = arctane 2 . 

Then a primitive function is 

F(x) = [ f(u) du = lmLi 2 (ie x ), 

J —oo 

(see [2H]) and the derivative is 

/'(*) = 



2 coshx 
Since 

S (P) = J2 {F{pk-Pj) + F{p3-Pk)-F{pj + Pk)-F{- Pj - Pk )--K{ Pj + Pk )) 



one obtains after some manipulations that 



= -2^(arctane Pfc Pl + arctan e Pk +Pl ) + $ r 



dpi 

J neighbors k 

This proves that Equations are the equations for a critical point of S(p). 

Equation (|ll)f) for the second derivative of S is obtained by taking the second 
derivative term by term in the first sum of equation l|14|l . For example, the 
second derivative of F(pk — Pj) is f'(pk — Pj){dpk — dpj) 2 . 

This concludes the proof of Proposition 01 
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